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ABSTRACT 

We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly 
discovered Milky Way satellite Reticulum IL Based on the spectra of 25 Ret II member stars selected 
from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 ± 0.5kms“® and 
a velocity dispersion of 3.3 ± 0.7kms“®. The mass-to-light ratio of Ret II within its half-light radius 
is 470 ± 2 IOM 0 /L 0 , demonstrating that it is a strongly dark matter-dominated system. Despite its 
spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC 
and SMC by 199 and 83 kms“®, respectively, suggesting that it is not gravitationally bound to the 
Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 
0.28 ±0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < —3. In combination 
with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. 

With a mean metallicity of [Fe/H] = —2.65 ± 0.07, Ret H matches Segue 1 as the most metal-poor 
galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight 
integral of the dark matter density squared is logjQ(J) = 18.8 ±0.6 GeV^ cm“^ within 0.2°, indicating 
that the predicted gamma-ray flux from dark matter annihilation in Ret H is lower than that of several 
other dwarf galaxies. 

Keywords: dark matter; galaxies: dwarf; galaxies: individual (Reticulum H); galaxies: stellar content; 

Local Group; stars: abundances 
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1. INTRODUCTION 


The population of known dwarf galaxies orbiting the 
Milky Way has grown rapidly over the last decade, with 
the discovery of the ultra-faint dwarfs by the Sloan Digi¬ 
tal Sky Survey (SDSS) more than d oubling the size of our 


Galaxy’s retinue of satellites (e.g., Will man et al. 
Zucker et al.|2006{||Belokurov et ah 2U(J7|. These ex' 


2005 


reme 


objects are the tocus of an enormous variety of on going 


work, ranging from their internal kinematics (e.g., Mar 

tin et al. 12007 

Simon & Geha|2007l 

Koposov et al. 

2011 

metallicities (e.g., Kirby et al.||2008 JNorris et al. 

2010 


Frebel et al.||2010| |Vargas et al.||2013 

to their star forma tion histories (e.g., 


h'rebel et al. 


Weisz et ai.l 


Brown et al. 20141, cosmological implications (e.g., Weisz 
et al.||2014b| r and ability to constrain dark matter mod¬ 
els viamdirect detection (e.g.,|Strigari et al.|2 007[|Acker- 
mannetal. 2014 Geringer-Saineth et al.|20M ). However, 
progress has slowed m the past hve years as th e flow of 
discoveries from SDSS dwindl ed ( although se e Lae vens 
et al.||2015 Martin et al.||2015 and |Kim et ai.||2015 ). 

'The re cent discovery ot eight new candidateaw arf 
galaxies (Bechtol et al. 2015 Koposov et al. 20151 in 
2000 deg^ of DEGam imaging d ata from the first year 


of the Dark Ener gy Survey (DES; Flaugher et al. 2015, 
Diehl et al.||2014) promises to reinvigorate studies ot the 


faintest galaxies. Some of the particularly interesting as¬ 
pects of these newly discovered objects are their apparent 
concentration around the Magellanic Glouds, the identi¬ 
fication of a relatively luminous dwarf (Eridanus II) near 
the Milky Way’s virial radius with possible recent star 
formation, and the existence of hyper-faint {My —3) 
dwarfs beyond the immediate vicinity of the Milky Way. 
Four of these new satellites — Reticulum II, Eridanus II, 
Tucana II, and Horologium I — can be fairly confidently 
classified as galaxies based on DES imaging alone, while 
the other systems have physical sizes and luminosities 
that overlap within their uncertainties with those of some 
globular clusters. 

In this paper we begin the process of spectroscopic 
follow-up observations of the new DES ultra-faint satel¬ 
lites. We present low and high resolution spectroscopy 
of stars in Ret II, the closest (d = 32 kpc) and best- 
characterized of the DES satellites, with an elliptic- 
ity o f 0.6 and a projecte d elliptical half-light radius of 
0.1° (Bechtol et al. 2015). The velocities and metallici- 
ties ot these stars conhrin that Ret II is indeed a dwarf 
galaxy. In ^we describe our spectroscopic target selec¬ 
tion, observations, and data reduction. We discuss the 


This paper includes data gathered with the 6.5 meter Mag¬ 
ellan Telescopes located at Las Campanas Observatory, Chile. 

f Based on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with the 
NSF on behalf of the Gemini partnership: the National Sci¬ 
ence Foundation (United States), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research Coun¬ 
cil (Australia), Ministerio da Ciencia, Tecnologia e Inovagao 
(Brazil) and Ministerio de Ciencia, Tecnologia e Innovacion Pro- 
ductiva (Argentina). 

t Based on data obtained from the ESO Science Archive Fa¬ 
cility under request number 157689. 


measurement of stellar velocities and metallicities and 
the classification of Ret II member stars in We focus 
on some of the implications of our results and examine 
how Ret II hts into the previously known population of 
Milky Way satellite galaxies in Q In ij^we summarize 
our hndings and conclude. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Target Selection 

Spectroscopic targets were selected from the object 
catalog derived from the coadded images of the first in¬ 
ternal annu al release of DES data (YlAl; Gruendl et 
al., in prep; ISevilla et al.||2011[ IDesai et al. |2012[ |Mohr| 


et al. 2012 Balbinot et al. 2015^. We identified ob¬ 
jectsrs stars based o n the spread model quantity out - 
put by SExtractor (Bertin 2011 Desai et al. 2012). 
Our stellar sample consists of well-rneasured objects with 
\spread_modelJ\ < 0.002 and flags-{g,r,i} < 4. We 
selected likely members based on a matched-filter maxi¬ 
mum likelihood procedure combining a spatial model of 
Ret II and the predicted color-magnitu de distribution 
of an old metal poor stellar population ([Bechtol et al. 


2015). We modeled the spatial distribution ot Ret II 


with the best-fit elliptical Plummer profile with ellip- 
ticity of 0.6 and an elliptical half-light radius of 0.1°. 
The distribution of stars in color-magnitude space was 
modeled by a composite of four isochrones bracketing a 


ties, Z = \ 

et al. 

2010|) 

ht dis 

tance 


-vr, 14.1 Gvrl, 

and metallici- 

([Marigo et al. 

[2008 

Girardif 


place: 


ance modulus of m — M = 17 . 5 and weighted by 


the initial mass function of Chabrier (2001). We used a 
composite isochrone to maximize the number of candi¬ 
date members while remaining agnostic to the exact age 
and metallicity of Ret II. 

One result of the maximum-likelihood fit is a mem- 
bers hip probability for each object in our stellar cata¬ 


log (Bechtol et al. 2015 Rozo et al. 2009). This mem- 
bersmp probability incorporates botn the spatial sepa¬ 
ration of the star from the centroid of Ret II and the 
distance of the object from the composite isochrone in 
color-magnitude space. Our spectroscopic follow-up sam¬ 
ple primarily consisted of stars with membership prob¬ 
ability p > 0.01 within 0.5° of the centroid of Ret II. 
Each of these objects was visually inspected in the DES 
coadded images and the imaging was confirmed to be of 
high quality. 

We targeted Ret II for follow-up spectroscopy with 
both Magellan/M2FS and Gemini/GMOS, and we have 
also made use of publicly available VLT/GIRAFFE spec¬ 
tra of Ret II stars in the ESO Archive. For observa tion s 
with the 256-fiber Magellan/M2FS spectrograph () |2.2| ), 
we targeted objects within 14.65' of Ret II with p > 0.01 
and g < 22. A small fraction of the 177 targets meeting 
these criteria were not observed because of fiber colli¬ 
sions and exclusion regions around bright stars. Given 
the availability of additional fibers, we also added 18 stars 
located near the Ret II isochrone that missed the proba¬ 
bility and/or photometric quality cuts. Out of these 195 
stars, a total of 185 were targeted with science fibers. Six 
fibers were broken, and the remaining 65 were placed on 
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blank sky positions chosen to have low count rates from 
a Digitized Sky Survey image. 


The Gemini/GMOS observation field (^2.3) was se¬ 
lected as the 5' x 5' region with the highest summed 
membership probability. Slits were iteratively priori¬ 
tized based on the membership probability and bright¬ 
ness of targets, focusing on stars on the giant branch for 
which metallicities could be measured via the Ga triplet 
(CaT) lines. Final mask creation was performed using 
the GMMPS v0.402 toolkit, which was used to assign 
slits to 33 science targets and 4 acquisitio n ob .iectsPj 

VLT/GIRAFFE observations of Ret II (j ]2.4^ were car- 
ried out under the auspices of the Gaia/ESO Survey 
(Gilmore et al. 20121, and the target selection was per- 
tormed by the Gaia/ESO Survey team independently of 
our DES photometric analysis chain. It is clear from the 
colors and magnitudes of the observed stars that they 
were chosen to lie near the Ret II fiducial sequence de¬ 
termined by Koposov et al. (2015), but we do not have 
access to the specihc criteria by which priorities were set. 
Comparing this sample with our catalog, we find that 14 
of the stars observed are high-probability member candi¬ 
dates {p > 0.5), 8 have low to intermediate membership 
probabilities (0.5 > p > 0.01), and 23 are unlikely to be 
members {p < 0.01). The majority (62) of the VLT spec¬ 
tra are of bright stars not included in the DES coadd cat¬ 
alog because of its sat uration limit at ^ 17, but (with 
one exception; (3.3.2) these stars appear to have been 
selected to be Milky Way disk stars and are not relevant 
for our purposes. 

2.2. Magellan/M2FS Spectra 

We obs erved Ret II with the multi-fiber M2FS spec¬ 
trograph (Mateo et al.|2012 1 on the Magellan/Clay Tele¬ 
scope on 2015 February 19. M2ES consists of two iden¬ 
tical spectrographs, generally referred to as “red” and 
“blue” even though neither is optimized for a particular 
wavelength range. Each spectrograph is fed by 128 1.2" 
fibers, which are positioned on the sky over a 29.5' diam¬ 
eter field with a plug plate. We used M2ES in its high 
resolution configuration at i? « 25000, with a narrow- 
band filter to isolate a single spectral order covering the 
Mg b spectral region (5120 — 5190 A) for each fiber. The 
256 M2FS fibers were placed on 185 stars selected from 
the DES YlAl internal data release, as well as 65 blan k 
sky positions for sky subtraction, as described in §2.1[ 

We obtained three 2400 s exposures on the Ret II field 
in variable seeing conditions and with decreased trans¬ 
parency resulting from thin clouds. We also obtained 
twilight sky spectra for flatfielding, ThArNe lamp frames 
for wavelength calibration, and spectra of the K giant ra¬ 
dial velocity standard star CD—43°2527 through a single 
fiber to check the velocity zero point. 

We reduced the data with the dohydra package in 
IRAF. We began by subtracting the bias level, recom¬ 
bining the files from the four amplifiers used to read out 
each CCD into a single frame, and masking out cosmi c 
rays using the L.A.Cosmic routine (van Dokkum 2001). 
We then trimmed the frames to the central ~ 1000 pix- 
els to eliminate contamination from neighboring spectral 
orders. We traced the fiber positions as a function of 
wavelength using a twilight frame and removed scattered 

http: //www. gemini . edu/node/12255?q=node/10458 


light with a two-dimensional fit to the areas between 
each group of 16 fibers. Because of the low signal-to- 
noise ratio (S/N) of the data (the brightest stars reach 
S/N « 20 pixel”and the large majority of the targets 
are at S/N < 10 pixel”^) and the excellent flatness of 
the e2V detectors (~ 1 %), we did not attempt to apply 
a flatfield correction (Massey & Hanson 2013). We ex¬ 
tracted the spectrum ot each hber in the 'I’hArNe lamp, 
twilight, standard star, and science frames. We also per¬ 
formed the same extraction on a variance image derived 
from the science frames to provide spectra of the pixel 
uncertainties. We determined the wavelength solution 
by fitting a fourth order polynomial to known Th and 
Ar line wavelengths on the ThArNe exposures. The po¬ 
sitions of the ThAr lines shifted by one pixel (3.8 kms”^) 
in the wavelength direction between the ThArNe frame 
taken immediately before Ret II was observed and the 
one taken immediately after. By examining the absorp¬ 
tion lines of the highest S/N stars, we determined that 
in the blue spectrograph, this shift occurred between the 
second and third Ret II exposures. In the red spectro¬ 
graph all three Ret II frames appeared to be aligned 
with each other. We therefore shifted the third frame 
taken with the blue spectrograph by one pixel and re-ran 
the reductions using only the ThArNe frame from before 
the Ret II observing sequence (to maintain consistency 
between the velocity scales for the red and blue spec¬ 
trographs, the wavelength solution for the red data was 
also determined with the ThArNe frame observed before 
Ret II). This wavelength solution was then applied to 
the data and the spectra were rebinned. We constructed 
master sky spectra for the red and blue spectrographs 
by coadding the ~ 30 sky fiber spectra obtained in each 
spectrograph. The sky spectra were scaled according to 
the throughput of each target fiber (determined from the 
twilight sky frames) and subtracted. 

2.3. Gemini/GMOS Spectra 

We a lso observed Ret I I with the GMOS-S spectro¬ 
graph (Hook et al. 2004) on the Gemini South Tele¬ 
scope in queue mode on several nights beginning on 2015 
February 15 through program GS-2014B-DD-8. We tar¬ 
geted 33 likely Ret H member stars with a multi-slit 
mask, using 0.75" slits and a minimum slit length of 
5". Our observations employed the R831 grating to pro¬ 
duce a spectral resolution of i? = 4400 over a range of 
2300 A per spectrum, and the RG610 filter to block sec¬ 
ond order light. A total of eleven 1200 s science-quality 
exposures were obtained. We observed at two different 
central wavelengths around the near-infrared CaT lines 
(8550 A and 8650 A) to ensure that key spectral features 
did not land in gaps between the three CCDs. Science 
observations were alternated with Quartz-Halogen flats 
and CuAr arc lamp exposures for calibration. 

We reduced the CMOS data using version 2.16 of the 
Gemini package in IRAF. For each dither on each night 
of observations, we created a flat field frame and used 
it to process the science frame(s). We combined expo¬ 
sures from the same dithers and then derived the wave¬ 
length solution. The spectra were rectified and trans¬ 
formed onto a common wavelength scale, and finally the 
dithers were coadded to eliminate the CCD gaps. We 
performed sky subtraction on each slit by masking out 
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the stellar continuum and fitting a linear function to the 
background in the spatial direction. Extractions were 
carried out on the sky-subtracted spectra. 


2.4. VLT/GIRAFFE Spectra 


Ret II was observed with the FLAMES/GIRAFFE 
spectrograph (Pasquini et al. 2000) on the VLT/UT2 


telesc ope as part ot the Gaia/ESU survey (Gilmore et al. 


20121. Because the Gaia/ESO survey is a public spectro- 


scopic survey, these data are public and available from 
the ESO Archiv^ immediately after the observations are 
obtained. 

Observations were taken in MEDUSA mode, which 
allows the simultaneous observation of up to 132 ob¬ 
jects. On the night of 2015 March 9, 2 x 1500 s expo¬ 
sure were taken using the HR21 grating, which covers 


the wavelength range from 8482-8981A at a resolution 
of i? ~ 16200. The same field was observed with a bluer 
wavelength setting in 2015 February, but because of the 
lower signal-to-noise ratio we do not make use of those 
data in this paper. The calibration frames taken as part 
of the Gaia/ESO Survey observations include 5 biases, 
3 flats, and 1 ThAr arc taken at the end of the night. We 
reduced the data with the GIRAFFE Gasgano pipeline 
(v2.4.8) provided by ESOQ This pipeline provides utili¬ 
ties for bias subtraction, flat-fielding, spectral extraction 
of individual objects and accurate wavelength calibra¬ 
tion. 

The GaT absorption lines are located close to a num¬ 
ber of bright sky emission lines, so we require accurate 
sky subtraction to recover their positions and equivalent 
widths. We first determined fiber to fiber wavelength 
offsets for the 16 GIRAFFE sky fibers using the night 
sky lines, finding that there are shifts of ^ 0.15 kms“^ 
between fibers. We then shifted the sky spectra based 
on these calculated offsets, linearly interpolated them, 
and combined them into a master sky spectrum, reject¬ 
ing outliers with a 3 (t clip. For each science spectrum, we 
determined the wavelength offset relative to the master 
sky and shifted the data accordingly. We scaled the sky 
spectrum to match the amplitude of the bright, isolated 
sky lines at 8885.85 A and 8919.6 A and subtracted it 
from the science spectrum. 


3. VELOCITY AND METALLICITY MEASUREMENTS 

Because the M2FS spectra have higher spectral resolu¬ 
tion than the other data sets and cover the largest sample 
of candidate Ret II member stars, they are particularly 
well suited for measuring the stellar kinematics of Ret II. 
The GIRAFFE spectra could also be used for accurat e 
radial velocity measurements (e.g., Koposov et al.|2011 ), 
but they include a smaller number of likely Ret 11 meni- 
bers, nearly all of which overlap with the M2FS target 
list. In order to avoid the complications of placing multi¬ 
ple instruments on a consistent velocity scale, we choose 
to use M2FS alone for the primary kinematic re sults re¬ 
ported in this paper. However, as we note in §4.1[ we 
obtain consistent results with the GIRAFFE velocities. 
The GMOS observations were made at much lower spec¬ 
tral resolution and would not be competitive with M2FS 


® http://archive.eso.org/cms/eso-data.html 
® http://www.eso.org/sci/software/gasgano.html 


or GIRAFFE for velocity work. GIRAFFE and GMOS 
velocities for stars not observed with M2FS are listed in 
Table but are not used to constrain the velocity dis¬ 
persion of Ret II. 

The wavelength range covered by the M2FS observa¬ 
tions is not well-suited to metallicity measurements us¬ 
ing the spectra alone, because the strength of the Mg b 
triplet lines depends on a number of stellar properties: ef¬ 
fective temperature, surface gravity, metallicity, and Mg 
abundance. No calibrations for the metallicities of indi¬ 
vidual stars as a function of Mg b strength are available 
in the literature. We therefore rely on the better-studied 
CaT region of the spectrum in the VLT and Gemini data 
for metallicity information. 

The primary data used in our analysis for the remain¬ 
der of the paper are 52 M2FS velocity measurements, 16 
GIRAFFE metallicities, and 6 GMOS metallicities. All 
velocity measurements presented here have been trans¬ 
formed to the heliocentric frame. We use heliocentric 
corrections of —8.26 kms“^ for the M2FS observations 
and —5.63 kms“^ for the GIRAFFE observations. 


3.1. Radial Velocity Measurements 

Given a velocity data set in which the intrinsic velocity 
dispersion may not be much larger than the uncertainties 
on individual velocity measurements, accurate character¬ 
ization of both the velocities and the measurement un¬ 
certainties is critical for a reliable determination of the 
velocity dispersion. We used Monte Carlo simulations 
to investigate the statistical uncertainties of our veloc¬ 
ity measurements, and repeat observations to constrain 
systematic uncertainties. 

We determi ned radial veloc i ties u sing the same basic 


approach as in Simon & Geha (2007) and subsequent pa¬ 
pers. We performed a fit to each spectrum with a high 


set between the two (|Lee et al. 12008 

Adelman-McCarthy 

et al. 2008 

Newman et al. 

2013). 

As a template spec- 


trum of the metal-poor giant HD 122563 obtained in 2009 
with Magellan/MIKE, because the M2FS RV standard 
CD—43° 2527 is relatively metal-rich and does not pro¬ 
vide a very good match to the spectra of the Ret H stars. 
We assume a heliocentric veloc ity of —26.5 kms“^ for 
HD 122563 (Chubak et al. 2012). We then ran 500 Monte 


Carlo simulations per star, m which noise comparable to 
the calculated variance in each pixel was added to the 
spectrum and the velocity was re-measured as above. 
The statistical uncertainty on the velocity of each star 
was defined to be the square root of the variance in the 
measured mean velocity for the 500 Monte Carlo spectra. 
Highly discrepant Monte Carlo iterations were discarded 
before computing the variance. The median statistical 
uncertainty on the velocity measurements for the M2FS 
data set is 1.0 kms“^. 

As a test of systematic velocity uncertainties in the 
M2FS data, we reduced the three Ret H science frames 
separately, and measured velocities for the brightest stars 
in each frame. Despite carrying out these reductions in 
multiple ways and computing velocities with two differ¬ 
ent techniques, we found that the velocity difference be¬ 
tween independent measurements of a star on frames i 
and j, Av = Vi — Vj, is consistently larger than would 
be expected from the statistical uncertainties on each in- 
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Figure 1. (a) DES color-magnitude diagram of Reticulum II. Stars within 14.65^ of the c enter of Ret II are plotted as small black dots, 
and stars selected for spectroscopy with M2FS, GIRAFFE, and CMOS (as described in EH are plotted as filled gray circles. Points 
surrounded by black outlines represent the stars for which we obtained successful velocity measurements, and those we identify as Ret II 
members are filled in with red. The four PARSEC isochrones used to determine membership probabilities are displayed as blac k lines, (b) 
Spatial distribution of the observed stars. Symbols are as in panel (a). The half-light radius of Ret II from|Bechtol et al.|(|2015|| is outlined 
as a black ellipse, (c) Radial velocity distribution of observed stars, combining all three spectroscopic data sets, i'ne clear narrow peak of 
stars at u ~ 60 kms“^ highlighted in red is the signature of Ret II. The hatched histogram indicates stars that are not members of Ret II; 
note that there are two bins containing non-member stars near u = 70 kms~^ that are over-plotted on top of the red histogram.. 


dividual velocity measurement {verr = + Verrj)- 

To obtain a standard deviation of Av/verr of one, we 
needed to add an additional 0.9 kms“^ uncertainty 
in quadrature with Verr- We consider this value of 
0.9 kms“^ to be the systematic uncertainty in the M2FS 
velocity measurements, and we define the total velocity 
uncertainty for each spectrum to be the quadrature sum 
of the Monte Carlo uncertainty from the template fit 
and the systematic uncertainty. A separate check on the 
systematics and their possible origin is provided by the 
high S/N M2FS twilight sky spectra. We fit these twi¬ 
light spectra with a high resolution solar template spec¬ 
trum. The scatter in velocity from hber to fiber was 
< 0.20 kms“^, so we conclude that the internal velocity 
errors over short timescales on an individual frame (in¬ 
corporating, e.g., any fiber-to-fiber systematics) are neg¬ 
ligible. However, over multiple science exposures span¬ 
ning several hours, this is not necessarily the case (see 
above). 

In order to verify the reliability of our velocity zero 
point, we also measured the velocity of the radial ve¬ 
locity standard star CD— 43°2527 by fitting it with the 
HD 122563 template, exactly as we did for the science 
spectra. For the two exposures on CD— 43°2527 , we find 
Vhei = 19.6 ± 0.1 kms“^ and Vhei = 19.9 ± 0.1 kms“^, 
compared t o the cataloged velocity of Vhei = 19.7 ± 
0.9 kms“^ ( Udry et'ar]|l999 ). 


3.2. Metallicity Measurements 
We calculated metallicities for 16 Ret H RGB stars 


with the CaT ca libration of Carrera et al. (2013). As 
recommended by Hendricks et af. "(2014), we measured 
the eq uivalent widths (FWs) of the CaT lines in the same 
way as |Carrera et al.[ fitting each of the three lin es with a 
Gauss i ^ plu s Lorentzian profile. Also following [Carrera] 
et al. ( 2013), we adopt the line and continuum regions 
defin ed by Cenarro et al. (2001), except for the 8498 A 


line. [Cenarro et af.| employed a continuum bandpass of 
8474 — 8484 A for this line, but the blue limit of the 


GIRAFFE spectra is 8482 A, so we instead use a region 
on the red side of the line from 8513 — 8522 A. This 
wavelength range may be modestly affected by two weak 
Fe I lines at 8514 A and 8515 A, but at the metallicity 
of typical ultra-faint dwarf stars any resulting depression 


of the continuum should be negligible over a 9 A band. 

CaT metallicity measurements usually use the horizon¬ 
tal branch (HB) magnitude to correct for the dependence 
of the CaT EWs on stellar luminosity. The horizontal 
branch magnitude of Ret H, however, is not well deter¬ 
mined because the galaxy contains so few HB stars. We 
therefore rely on the calibratio n of CaT EW a s a fun ction 
of absolute V magnitude from Carrera et al. (2013). We 
convert the DES g and r magnitudes to the SDSS photo- 
metric syst em, and then use th e relations for metal-poor 
stars from Jordi et al. (2006) to transform to V. We 
determine absolute magnitude s assuming a distance of 
32 ± 3 kpc ( Bechtol et al. 12015 ) and a V-band e xtinction 
of Ay = 0.05 mag ( [Schlaify fe h'inkbeiner||2011 ). 


3.3. Spectroscopic Membership Determination 
3.3.1. M2FS 


Out of the 185 M2FS fibers placed on stars, we success¬ 
fully measured velocities for 52, including a large major¬ 
ity of the observed targets brighter than g = 20.6. The 
remaining stars had S/N ratios too low for spectral fea¬ 
tures to be confidently detected in the data. The veloc¬ 
ity measurements and other properties of the stars are 
listed in Table The velocity distribution we measure 
from the M2FS spectra exhibits a strong peak at a veloc¬ 
ity of ~60 kms“^ (see Fig. [^, as is characteristic of a 
gravitationally bound system. Approximately half of the 
stars for which we measure velocities are contained in this 
peak, with the remainder spread across a wide range of 
heliocentric velocities from ~0 kms“^ to ~330 kms“^. 

For a large majority of the observed stars, the member¬ 
ship status is unambiguous; stars with Vhei > 90 kms“^ 
and Vhei < 40 kms“^ are clearly not related to the peak 
associated with Ret H, while those very near the mean 
velocity of the system and close to the central position 
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spatially are almost certainly members. However, to en¬ 
sure that the member sample is defined optimally we 
carefully examine all stars within 20 kms“^ of the mean 
velocity of Ret II, considering their velocities, positions in 
the color-magnitude di agram, spatial locati ons, member¬ 
ship probabilities from Bechtol et al. (2015), and spectral 
features. Below we discuss the individual stars whose 
membership is not immediately obvious. 

Three stars in our sample have velocities of 
kms“^, just to the left of the Ret II peak in 
Fig. [^, and about 15 kms ^ away from the systemic 
velocity. Of these three, DES J033405.49—540349.9 
is very metal-rich, DES J033437.34—535354.0 is a 
double-lined spectroscopic binary with two strong cross¬ 
correlation peaks separated by ^ 60 kms“^, and 

DES J033540.70—541005.1 is located close to two half- 
light radii away from the center of Ret II along the minor 
axis and is separated from the mean velocity by ~4(T. 
None of these stars has characteristics consistent with 
membership in Ret II. DES J033540.70—541005.1 has 
colors consistent with what would be expected for a red 
horizontal branch star in Ret II, but is ~ 0.1 mag fainter 
than the isochrone shown in Eig. Given its spatial 
and velocity offsets from Ret II as well as its position in 
the color-magnitude diagram, it is almost certainly not 
a member. Spectra of these three non-member stars are 
displayed in Fig. 

Another group of four stars is located slightly to 
the blue of the isochrones we use to describe the 
Ret II stellar population (by ^ 0.02 — 0.06 mag). 
DES J033524.00-540226.7 is the farthest from the 


isochrone, but is also the faintest of these stars (sitting 
near the base of the red giant branch) and consequently 
has the largest photometric uncertainties. Its velocity is 
within 3 kms“^ of the systemic velocity of Ret II and it 
is located inside the half-light radius, so we consider it 
a member. DES J033635.78—540120.2 is also relatively 
faint, but is farther away from Ret II in both position (at 
^ 1.3i?i/2) and velocity (6 kms“^). It is likely a member 
of Ret II, but the spatial and velocity offsets make that 
classification less certain. DES J033531.14—540148.3 is 
one of the brightest candidate members, in the clump of 
stars near g ~ 17.5 in Eig. Ba- It is within the half-light 
radius of Ret II and has a very low metallicity as deter¬ 
mined from the GIRAFFE GaT measurements. Because 
of its low metallicity and spatial position, we conclude 
that it is a Ret B member. The color offset to the blue 
of the giant branch could indicate that it is an asymp¬ 
totic giant branch star. DES J033550.10—540139.2 is 
located quite close to the center of Ret B, but is one of 
the bluest (non-HB) stars with a velocity consistent with 
that of Ret B. This star is projected quite close to the 
center of Ret B and has extremely weak absorption lines. 
If it is indeed a member, it could have an unusually low 
metallicity. As a result of the weak spectral features, the 
velocity is quite uncertain; it is consistent with member¬ 
ship, but it has little impact on the derived properties of 
Ret B because of the large uncertainty. 

Two stars are offset to the red side of the isochrone. 
DES J033436.70—540645.0 is almost two half-light radii 
from the center of Ret B, although it is very close to the 
position of the velocity peak. It also has strong absorp¬ 
tion lines, indicating a relatively high metallicity. Given 



Figure 2. Magellan/M2FS spectra in the Mg b triplet region for 
the three stars near the edge of the Ret II velocity distribution. 
The wavelengths of two Mg lines and an Fe line are marked in the 
bottom panel, and the third component of the Mg triplet is just 
visible at a wavelength of 5185 A at the right edge of each spec¬ 
trum. The spectrum of DES J033540.70—541005.1 (top) appears 
similar to that of a Ret II member, but the color, spatial position, 
and velocity offset of this star make that classification unlikely. 
The very strong Mg absorption in DES J033405.49—540349.9 (mid¬ 
dle), as well as the wealth of other absorption features on the 
blue side of the spectrum, indicate that the star is more metal- 
rich than would be expected for a system as small as Ret II. 
DES J033437.34—535354.0 (bottom) is a double-lined binary star 
with a velocity separation of ~ 60 kms“^. The redshifted ab¬ 
sorption component from the secondary star is most visible in the 
middle line of the Mg triplet. 

that it is offset from the rest of the Ret B population in 
color, spatial position, and metallicity, we conclude that 
it is probably not a member. DES J033544.18—540150.0 
is at the centroid of the Ret B velocity peak within the 
uncertainties, and is also located very close to the central 
position, but is ~ 0.04 mag redder than the isochrone. 
Since it is near the base of the giant branch, the pho¬ 
tometric uncertainties could contribute to this offset in 
color, and we consider DES J033544.18—540150.0 a likely 
member of Ret B. 

Because the stars for which membership is plausible 
have velocities quite similar to that of Ret B (and in 
some cases have large uncertainties), including or exclud¬ 
ing them from the member sample does not have any 
significant effect on the properties we derive for Ret B in 
^ We show the correspondence between M2FS spectro¬ 
scopic members and high photometric membership prob¬ 
ability in Eig. 1^ 

3.3.2. GIRAFFE and CMOS 

We also identify a handful of Ret B members in the 
GIRAFFE and GMOS data sets that were not observed 
with M2ES. We use a velocity measurement based on the 
Paschen lines to confirm tha t the candidate BHB star 
DES J033539.85-540458.1 (^^ observed by GMOS 
is indeed a member of Ret B, with a velocity of 69 ± 
6 kms“^. The GIRAFEE targets included a bright {g ~ 
16.5) star at ( 02000 :^ 2000 ) = (03:35:23.85,-54:04:07.5) 
that was omitted from our photometric catalog and 
M2FS observations because it is saturated in the coad- 
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ded DES images. However, the spectrum of the star 
makes clear that it is very metal-poor and is within a 
few kms“^ of the systemic velocity of Ret II. While the 
magnitudes derived from individual DES frames place it 
slightly redder than the isochrone that best matches the 
lower red giant branch of Ret II, it is also located inside 
the half-light radius, and is very likely a member. In fact, 
it is probably the brightest star in any of the ultra-faint 
dwarfs. 

Another GIRAEFE target, DES J033548.04-540349.8, 
was not included in the M2FS sample because it was as¬ 
signed a photometric membership probability of zero by 
the maximum-likelihood analysis (|Bechtol et al. 2015). 
It is very slightly redder (relative to the isochrone) than 
the bulk of the Ret II stars, but it is located within the 
half-light radius and has a velocity ^ 1 kms“^ from the 
systemic velocity. The metallicity we determine from the 
GIRAEFE spectrum is low, but somewhat more metal- 
rich than any of the other member stars, perhaps con- 
sis tent with its color. The Mg I A8807 line advocated 
by|Battaglia & Starkenburg| (|2012|) as a discriminant be- 
tween dwarf-galaxy red giants and foreground main se¬ 
quence stars is weak (EW = 0.1 A), suggesting that it 
is a member of Ret B, but since it would be the highest 
metallicity star in the galaxy we regard this conclusion 
as tentative. Although all three of these stars are likely 
members of Ret B, we do not include them in the de¬ 
termination of the kinematic properties in §4.1| because 
they have not been placed accurately enough in the same 
velocity reference frame. 

We classify the GIRAFFE target 
DES J033524.96—540230.7 as a non-member based 
on its extremely broad CaT lines, despite its overlap 
with the Ret B velocity distribution. 

Finally, we note that for DES J033515.17—540843.0, 
the M2FS and GIRAFFE velocity measurements differ 
by 8.2±3.5 kms“^. The significance of this difference is 
2.4cr, corresponding to a p-value of 2%. While we cannot 
rule out the possibility that this difference in velocity is 
attributable to chance or systematics, we suggest that 
this star may be a binary system in Ret B. Orbital mo¬ 
tion of 8 kms“^ over a time span of 18 days is consistent 
with the handful of pr eviously detected binaries in ultra- 
faint dwarf galaxies (Koposov et al. 2011 Koch et al. 
20T4|. 


3.3.3. Foreground Modeling 

Finally, we consider the expected contamination from 


sancon (Robin et al. 

2003D and TRILEGAL (Girardi 

et al. 

20U5 |Vanholle 

beke et al. 

2009) Galactic stellar 


foreground stars in the region of Ret B. We se lect stars 


within 0.02 mag of the composite isochrone ((2.1) and 
with n < g < 20.6. From the resulting velocity distri- 
bution, we find that ~ 10% of foreground stars passing 
these photometric cuts have a velocity consistent with 
the heliocentric velocity peak of Ret B (55 - 70 kms“^). 
We set the overall normalization for the number of fore¬ 
ground stars by assuming that all stars in the sideband 
of the peak {v < 40kms“^ or u > 90kms“^) belong 
to the foreground distribution. The expected number of 
foreground stars within the peak of ranges from 3.0 (Be- 
sangon) to 4.1 (TRILEGAL). Given the small number 



Figure 3. Comparison of photometric membership probabilities 
determined from a maximum-likelihood fit to the DES data and 
spectroscopic membership as determined from the velocity mea¬ 
sured by M2FS. 

statistics, this is consistent with the results of the mem¬ 
bership analysis above. Under these assumptions, the 
foreground-only hypothesis for the Ret B velocity peak 
is excluded at p « 10“^^. 

3.4. Horizontal Branch Stars 

The DES photometric catalog contains three can¬ 
didate blue HB (BHB) stars, three candidate red 
HB stars, and one star located slightly below the 
HB at an intermediate color that could be consis¬ 
tent with the RR Lyrae instability strip. We tar¬ 
geted all seven of these stars with M2FS, but were 
only able to determine velocities for five of them. For 
the BHB stars, we confirm DES J033618.68-535745.1 
as a member of Ret B, but we do not detect any 
spectral features for DES J033539.85—540458.1, and 
DES J033612.7—535602.2 was observed with a fiber 
with 3% throughput so the spectrum was not use¬ 
ful. As mentioned in §3.3[ our GMOS spectrum of 
DES J033539.85—540458.1 confirms it as a Ret B mem¬ 
ber. The velocities of the three RHB candidates indicate 
that they are not members. The final HB candidate, 
DES J033341.71—540007.3, has a measured velocity of 
34.4± 1.2 kms“^, which could be consistent with mem¬ 
bership in Ret B if it were an RR Lyrae variable, but 
we do not detect photometric variations for it or any of 
the BHB stars in the 2 — 4 DES single-epoch images per 
filter available in this field. 


4. DISCUSSION 

4.1. Velocity Dispersion and Mass 

After identifying likely members of Ret B with a combi¬ 
nation of photometry, velocity, and metallicity informa¬ 
tion, we can calculate the systemic velocity and velocity 
dispersion of Ret B. For the sample of 25 member stars 
with M2FS velocity measurements identified in ^ we 
find a systemic vel ocity of Vhei = 62.8 ± 0.5kms“ . Us¬ 
ing the method of Walker et al. (2006), we determine a 
maximum-likelihood estimate tor the velocity dispersion 
of cr„ = 3.3 ± 0.7kms“^. If we instead rely on the GI¬ 
RAFFE data set (which includes 17 member stars), we 
obtain fully compatible values of Vhei = 63.3 ±0.8 km s“^ 
and (T„ = 3.5 ± 1.0kms“^. The velocity dispersion of 
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Figure 4. Mass-to-light ratio (left) and metallicity (right) of Ret II compared to other galaxi es within 1 Mpc. Ma ss-to-light ratios were 
calculated from the velocity dispersions, half-light radii, and absolute magnitudes compil ed bvlMcConnachiel (I2012t in order to maintain 
consi stency with t he va lue we report for Ret II. Metallicities were also adopted from the |McConnachj^ l|iiUi2|| compilation, with updates 
from [Kirby et aL]||2013|l where newer measurements were available. 


Ret II is also robust to membership uncertainties; deci¬ 
sions about includi ng o r excluding the plausible member 
stars discussed in §3.3| change the calculated dispersion 
by less than the stated uncertainty. Only by adding stars 
that are very unlikely to be members for multiple rea¬ 
sons can the best-fit velocity dispersion be significantly 
increased. In principle, the low velocity dispersion of 
Ret II increases the likelihood that the measured disper¬ 
sion could be artificially inflated because of the orbital 
motions of binary stars. While our observations do not 
span a long enough time baseline to have much sensitiv¬ 
ity to binaries, previous studies have shown that even in 
the smallest dwarfs binary stars do n ot significantly in- 


Martinez et al. 

2011 

Wolf et al. (P 

010| 


locity anisotropy of a stellar system, the mass enclosed 
within the half-light radius can be accurately computed 
as 

^2 \ /R, 


A^x/2 — 930 


km2 s 2 


fl/2 

pc 


Mq. 


( 1 ) 


Using this relation, the mass within the elliptical half- 
light radius of Ret II {R 1/2 = 55 ± 5 pc; Bechtol et al. 
20151 is 5.6 ± 2.4 x 10 ^ Mf^. The ab s olute magnitude of 
My = —3.6 ±0.1 from [Bechtol et al. (20151 translates to 
a luminosity of 2360 L©, leading to a mass-to-light ratio 
of 470 ± 2 IOM 0 /L 0 . This value is consistent with the 
inverse correlation between mass-to-light ratio and lumi¬ 
nosity for other Local Group dwarf galaxies, although 
Ret II is on the low end for galaxies of similar luminosity 
(see Fig. Ill . 

4.2. Metallicity Spread 

Even a visual inspection of the M2FS spectra reveals 
clear chemical abundance differences within the Ret II 
member sample. In Fig. [5] we display spectra of three of 
the stars within the clump of Ret II members at g « 18.9 
in Fig. The absolute strength of the Mg absorption 


and the Fe absorption, as well as the ratio between the 
two, varies from star to star, demonstrating that Ret II 
is not chemically homogeneous. 

The CaT metallicity measurements from the 
VLT/GIRAFFE spectra provide a more quantita¬ 
tive demonstration of the spread in overall metallicity 
from star to star in Ret II. We find member stars 
spanning more than 1 dex in iron abundance, ranging 
from [Fe/H] = —2.0 to [Fe/H] = —3.3. The individual 
metallicities are displayed as a function of radius in 
Fig. i!^. There is no indication of a metallicity gradient 
within Ret II. The metallicity spread, calculated with 
the same technique we used in §4.1| for the velocity 
dispersion, is 0. 28 ± 0.09 dex. It we instead employ 


the definition of Kirby et al. (2011), we find a similar 
spread. Both this dispersion and the presence of several 
extremely metal-poor stars at [Fe/H] < —3 demonstrate 
that Ret II has the chemical characteristics of a dwarf 
galaxy rather than a globular cluster. The mean 
metallicity of Ret II is [Fe/H] = —2.65 ± 0.07, ranking 
it along with Segue 1 as the most metal-poor galaxies 
known. This metallicity measurement places Ret H 
somewhat belo w the stellar mass-m etallicity relationship 
determined by [Kirby et al.j (|2013|) and argues against 
the existence of a metallicity floor for dwarf galaxies 
(Fig. [4|. The kinematic and chemical properties of 
Ret Ilare summarized in Tabled 

4.3. Reticulum II and the Magellanic Clouds 

The proximity of the DES satellites to the Magel¬ 
lanic Clouds suggests that some or all of them may have 
originated in a Magellanic group. Ret H is the clos¬ 
est of these objects to the Magellanic Clouds, with a 
three-dimensional separation of just 23 ± 3 kpc. The ra¬ 
dial velocity of Ret H differs from that of the LMC by 
200 kms“^, which represents a lower limit on the dif¬ 
ference in their space motions. The total mass of the 
LMC out to a distance of 23 kpc is not known, but 
Kallivayalil et al. (2013) argue that reasonable bounds 
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Figure 5. (Left) Magellan/M2FS spectra in the Mg b triplet region for three Ret II member stars covering a range of line strengths. From 
top to bottom, the stars are DES J033556.28—540316.3, DES J033454.24—540558.0, and DES J033457.57—540531.4. These stars span 
only 0.1 mag in luminosity and 0.08 mag in g — r color, so their effective temperatures and surface gravities should be very similar. Any 
differences in line strength therefore translate directly into chemical abundance differ ences. The apparen t emission features near 5182 A in 
the spectrum of DES J033454.24—540558.0 are contamination by the Littrow ghost (|Burgh et al.|2007||. (Right) VLT/GIRAFFE spectra 
of the bluer two CaT lines for the same stars. 


Table 1 

Summary of Properties of Reticulum II 


Row 

Quantity 

Value 

(1) 

RA (J2000) 

03:35:41 

(2) 

Dec (J2000) 

-54:03:00 

(3) 

Distance (kpc) 

32 

(4) 

LIvfi 

-3.6 ±0.1 

(5) 

Lv,o (Lq) 

2360 ± 200 

(6) 

€ 

U.DU_q 20 

(7) 

ri/2 (pc) 

55 ±5 

(8) 

Vhel (km8“^) 

62.8 ± 0.5 

(9) 

Vgsr (kms“U 

-92.5 ±0.5 

(10) 

a (kms~^) 

3.3 ±0.7 

(11) 

Mass (Mq) 

5.6 ±2.4 X IQS 

(12) 

M/Lv (M 0 /L 0 ) 

470± 210 

(13) 

Mean [Fe/H] 

-2.65 ± 0.07 

(14) 

Metallicity dispersion (dex) 

0.28 ±0.09 

(15) 

login -1(0.2°) (GeV2 cm-5) 

18.8 ±0.6 

(16) 

login •7(0.5°) (GeV^ cm-®) 

18.9 ±0.6 


Note. — Rows (l )-(7) are taken f rom th e DES pho¬ 
tometric analysis of|Bechtol et al.| ||2015||. Values in 
rows (8)-(16) are derived m this paper. 


are 3 x 10^*^ M,^ to 2 5 x 10^^ Mq. Using DES imaging, 


Balbinot et al. (20151 determine a truncation radius of 
13.5 ± 0.08 kpc for the LMC. If this truncation is tidal in 
origin, the fa vored mass is at the lower limit of the range 
suggested by Kallivayalil et al. (2013). For a mass of 
3x 10^° Mq, the escape velocity of the LMC at d = 23 kpc 
is only 106 kms“^, and for Mlmc < 10^^ Mq the escape 
velocity is less than 200kms“^. Only if the LMC mass 
is higher than this value and the tangential velocity of 
Ret II is similar to that of the LMC could Ret II cur¬ 
rently be a bound satellite of the Magellanic Clouds. Of 
course, even if Ret II is not bound to the LMC now, it 


is possible that it was an LMC satellite at earlier times 
and was stripped from that host by the tidal field of the 
Milky Way more recently. A proper motion measure¬ 
ment of Ret II will be needed to conclusively determine 
its origin, and its location just 32 kpc away from the Sun 
should enable the proper motion to be determined in the 
next few years. 

The observed velocity distribution of high velocity 
{vhei > 100 kms“^) stars displayed in Figure ap¬ 
pears inconsistent with the distribution of veloci ties fo r 
Milky Way stars predicted by Galactic models (p.3.3). 
In particular, there are unexpected concentrations of 
stars around 145 kms“^ (six stars) and 300 kms“^ (five 
stars). The former group of stars is near the heliocen¬ 
tric v elocity of the SMC (145.6kms“^; Harris & Zaritsky 


2006|) , suggesting that these might be distant SMC stars. 


It they are at the distance of the SMC, they are at pro¬ 
jected separations of 27 kpc, approximately twice as far 
out as the most distant currently known SMC popula¬ 
tion ( Nidever et al.|20fT ), and have therefore likely been 
tidally stripped, 'i’he higher velocity stars have very sim¬ 
ilar velocities to the Magellanic Stream gas a few degrees 
away from Ret II, and could represent a stellar counter¬ 
part of the Stream. 


4.4. J-Factor 

It is posited that dark matter par ticles could self- 


annihilate to produce gamma rays (e.g., Gunn et a l. 1978. 
Bergstrom & Snellman 1988 iialtz et al. 2U(J8[ ).' 'The 


large dark matter content, relative proxirnity, and low 
astrophysical foregrounds of dwarf galaxies make them 
promising targets for the detection of these gamma rays. 
The predicted signal from the annihilation of dark mat¬ 
ter particles is proportional to the line-of-sight in tegral 


through the square of the dark matter density (e.g., Baltz 
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eFaI][200^ , 


J(Arj)= f [ 

J Af2 J 1.0.s. 


( 2 ) 


Here, /9 dm(»’) is the dark matter particle density, and the 
integral is performed over a solid angle AH. The J-factor 
is derived by modeling the velocities using the spheri¬ 
cal Jeans equation, with assumptions on the theoretical 


halo (e.g.. 

Striga 

ri et al.|2008 Essig et al.|2009 Charbon- 

nier et al. 

|20II[ |Martinez||20I3| [Geringer-Sameth et al. 

20I5ap. Here, we model the dark matter halo as a gen- 


et al. 1997). We use flat, ‘uninfo rmative’ priors on the 


dark matter halo parameters (see Essig et al.||2009) and 


assume a constant stellar velocity anisotropy. Using this 
procedure, we find an integrated J-factor for Ret 11 of 
logio('^) = ^ 0-6 GeV^ cm“® within an angular cone 

of radius 0.2°, and logio(J) = 18.9 ± 0.6 GeV^ cm“^ 
within 0.5°. This latter value assumes that the dark mat¬ 
ter halo extends beyond the radius of the outermost spec¬ 
troscopically confirmed star, but truncates within the es¬ 
timated tidal radius for the dark matter halo 1 kpc). 
The quoted uncertainties are Icr, and are estimated by 
modeling the posterior probability density function of 
logio(<J) ^ Gaussian. Note that the uncertainty ob¬ 
tained by modeling this individual system is larger than 
is obtained by m odeling the entire population of dSphs 
( Martinez||2013 ). 

Several previously known ultra-faint dwarf galaxies 
possess larger mean J-factors than Ret 11, most notably 
Segue 1, Ursa Major II, and Goma Berenices ([Acker- ' 
maim et al.||2014| [Geringer-Sameth et al.||2015a||(Jonrad 
et al.||2015p. Though the velocity dispersions of Ret 11 
and Segue 1 are consistent within uncertainties. Ret II 
is more distant (32 kpc compared to 23 kpc) and has a 
larger half-light radius as measured along the major axis 
(55 pc compared to 29 pc). The larger distance and larger 
half-light radius imply a reduced mean J-factor relative 
to Segue 1. In comparison to Ursa Major II, Ret II is 
at a similar distance, but has a velocity dispersion that 
is smaller by roughly a factor of two. The larger dis¬ 
persion, and hence mass, accounts for the larger J-factor 
of Ursa Major II. Coma Berenices is more distant than 
Ret B (44kpc compared to 32kpc); however, the larger 
velocity dispersion of Coma Berenices implies a slightly 
larger mean J-factor. 

Since Segue 1, Ursa Major B, and Coma Berenices 
all possess larger J-factors than Ret B, we expect dark 
matter annihilation to produce a larger gamma-ray flux 
from these objects. However, no gamma-ray excess has 
been asso ciated with any of the pr eviously known dwarf 
galaxies (Ackermann et al. 2015). Given comparable 
gamma-ray sensitivity, it is unlikely that a dark matter 
annihilation signal would be detected from Ret B with¬ 
out also b eing detected from dwarf ga l axies with higher 
J-factors dPrlica- Wagner et al . [201 5 Geringer-Sameth 
et al.||20I5D ' Hooper fc Linden||2UI5[ ). 


5. SUMMARY AND CONCLUSIONS 

We have presented the first spectroscopic analysis of 
the recently discovered Milky Way satellite Reticulum B. 
We measure the velocities of 25 Ret B members us¬ 


ing high resolution spectroscopy from Magellan/M2FS, 
as well as metallicities from the Ca triplet lines for 
16 members using high resolution spectroscopy from 
VLT/GIRAFFE and 6 members using low resolution 
spectroscopy from Gemini South/GMOS. Ret B has 
a velocity dispersion of (t„ = 3.3 ± 0.7kms“^, corre¬ 
sponding to a dynamical mass within its half-light ra¬ 
dius of 5.6 ± 2.4 X 10^ Mq and a mass-to-light ratio of 


47O±2IOM0/L0. 

The metallicity of Ret B determined from the GaT is 
[Fe/H] = —2.65 ± 0fo 7, consistent with that of Segue I 


2014) within the uncertainties, and 


(Frebel et al. 

0.2 dex lower than that of any other known galaxy. We 
find that Ret B has an internal metallicity spread of 
0.28±0.09 dex, with stars spanning a total range of more 
than I dex. Even the most metal-rich stars in the galaxy 
are at [Fe/H] « —2.0, and Ret B contains 4 extremely 
metal-poor stars with [Fe/H] < —3. Both its chemical 
and kinematic properties confirm that Ret B is a dwarf 
galaxy. 

The location of Ret B just 23 kpc away from the LMC 
suggests that it could have originated as a satellite of 
the Magellanic system rather t han having always been 
associated with the Milky Way (Bechtol et al.||20I5 Ko- 
posov et al. 2015). However, our measured systemic 
radial velocity of Vhei = 62.8 ± 0.5kms“^ means that 
Ret B is moving away from the LMC at a minimum ve¬ 
locity of I99kms“^. According to current LMC mass 
estimates, this velocity likely exceeds the escape veloc¬ 
ity of the LMC, indicating that the two objects are not 
gravitationally bound. This result does not rule out the 
possibility that Ret B was previously a Magellanic satel¬ 
lite, and future proper motion measurements will shed 
more light on its origin. 

The J-factor calculated from the internal kinematics of 
Ret B is logio ^ =*= 0.6 GeV^ cm“® within a radius 

of 0.2°, somewhat lower than prev iously estimated based 
on th e galaxy’s distance alone (Drlica-Wagner et al. 
2015). The predicted gamma-ray hux from dark mat- 


ter annihilation in Ret B is therefore likely to be lower 
than that predicted for several other Milky Way satel¬ 
lites ([Dr lica- Wagner et ah ] 2015 Geringer-Sameth et al. 


2015b} ' Hooper &: Linden zOISj ). 
Satellite galaxies like Ret II 


galaxies like Ret B provide a crucial testing 
ground for the ACDM paradigm, and accordingly, the 
search for ultra-faint galaxies has become a major theme 
of near-held cosmology. It is expected that many ad¬ 
ditional Milky Way satellite galaxies could be found in 
ongoi ng and near-future wide-held optical im aging sur¬ 


veys (Tollerud et al. 2008} Hargis et al. 2014). The link 
between newly discovered stellar systems and the dark 
matter halos in which they may reside is established by 
follow-up dynamical and chemical analysis. Ret B is the 


et al.|20I5l|Ko] 

□osov et al.|20I5MLaevens et al.|20I5||Mar-| 

tin et al. [2015 

Kim et al. 2015) to be spectroscopically 


lite galaxy. The spectroscopic campaign to characterize 
new satellite galaxy candidates represents an essential 
step in ongoing tests of the standard cosmological model. 


This publication is based upon work supported by the 
National Science Foundation under grant AST-1108811. 
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Table 2 

Velocity and metallicity measurements. 


ID 

ra 

DEC 

9 ^ 

V 

EWgiraffe [Fe/HjoiRAFFE 

EWgmos 

[Fe/HjoMOS MEM 


(deg) 

(deg) 

(mag) (mag) 

( kms“^) 





M2FS 


DES J033341.71-540007.3 

53.42377 

-54.00201 

18.26 

18.20 

34.42 ± 1.16 

- 

- 

- 

- 

NM 

DES J033405.49-540349.9 

53.52287 

-54.06387 

17.58 

16.99 

50.96 ± 0.94 

- 

- 

- 

- 

NM 

DES J033413.06-535956.1 

53.55441 

-53.99892 

19.45 

18.99 

196.81 ± 1.44 

- 

- 

- 

- 

NM 

DES J033413.94-540934.4 

53.55807 

-54.15956 

18.26 

17.68 

16.70 ± 1.06 

- 

- 

- 

- 

NM 

DES J033418.32-541006.2 

53.57632 

-54.16839 

18.02 

17.45 

24.41 ± 0.98 

- 

- 

- 

- 

NM 

DES J033429.94-541111.8 

53.62474 

-54.18661 

19.77 

19.32 

286.61 ± 1.55 

- 

- 

- 

- 

NM 

DES J033430.07-540922.2 

53.62528 

-54.15618 

18.17 

17.85 

267.19 ± 0.97 

- 

- 

- 

- 

NM 

DES J033436.70-540645.0 

53.65291 

-54.11249 

18.85 

18.29 

70.52 ± 1.01 

- 

- 

- 

- 

NM? 

DES J033437.34-535354.0 

53.65557 

-53.89832 

17.69 

17.11 

52.94 ± 1.10 

- 

- 

- 

- 

NM 

DES J033437.98-541359.7 

53.65825 

-54.23324 

20.55 

20.07 

137.64 ± 2.58 

- 

- 

- 

- 

NM 

DES J033439.66-540754.4 

53.66523 

-54.13177 

19.38 

18.87 

59.32 ± 1.84 

- 

- 

- 

- 

M 

DES J033439.81-540058.5 

53.66587 

-54.01626 

18.21 

17.93 

197.28 ± 1.00 

- 

- 

- 

- 

NM 

DES J033447.94-540525.0 

53.69974 

-54.09028 

17.52 

16.92 

62.28 ± 1.03 

- 

- 

- 

- 

M 

DES J033449.20-535019.7 

53.70502 

-53.83881 

17.62 

17.02 

13.13 ± 0.93 

- 

- 

- 

- 

NM 

DES J033453.23-541403.9 

53.72180 

-54.23443 

20.59 

20.19 

285.07 ±2.73 

- 

- 

- 

- 

NM 

DES J033453.50-540454.6 

53.72291 

-54.08182 

20.18 

19.80 

133.74 ± 2.58 

- 

- 

- 

- 

NM 

DES J033454.24-540558.0 

53.72600 

-54.09945 

18.95 

18.44 

69.70 ± 1.38 

1.69 ± 0.25 

-2.78± 0.17 

- 

- 

M 

DES J033457.57-540531.4 

53.73988 

-54.09206 

18.94 

18.42 

59.84 ± 1.17 

2.77 ± 0.23 

-2.19 ± 0.12 

- 

- 

M 

DES J033502.50-540354.3 

53.76041 

-54.06507 

19.24 

18.75 

67.94 ± 1.14 

1.98 ± 0.43 

-2.54 ± 0.26 

- 

- 

M 

DES J033502.87-540109.8 

53.76196 

-54.01940 

20.36 

19.93 

67.96 ± 3.48 

- 

- 

- 

- 

M 

DES J033506.56-540604.3 

53.77734 

-54.10120 

19.79 

19.28 

158.97 ± 1.36 

- 

- 

- 

- 

NM 

DES J033509.50-540229.7 

53.78959 

-54.04158 

18.24 

17.76 

95.55 ± 0.96 

5.72 ± 0.29 

- 

- 

- 

NM 

DES J033511.66-540321.8 

53.79858 

-54.05606 

19.31 

18.83 

65.64 ± 1.34 

1.27± 0.27 

-3.02 ± 0.25 

- 

- 

M 

DES J033513.73-540456.7 

53.80722 

-54.08242 

19.70 

19.24 

57.43 ± 2.38 

2.26 ± 0.40 

-2.29 ± 0.23 

- 

- 

M 

DES J033514.01-540558.2 

53.80839 

-54.09949 

20.01 

19.57 

63.47 ± 1.45 

- 

- 

- 

- 

M 

DES J033515.17-540843.0 

53.81322 

-54.14529 

19.73 

19.25 

67.87 ± 1.36 

2.19± 0.58 

-2.33 ± 0.33 

- 

- 

M 

DES J033517.01-540403.0 

53.82088 

-54.06751 

19.72 

19.24 

66.31 ± 1.40 

0.93 ± 0.38 

-3.31 ± 0.54 

1.37 ±0.20 

-2.86 ± 0.18 

M 

DES J033520.37-541816.9 

53.83486 

-54.30470 

20.53 

20.18 

311.02 ± 2.83 

- 

- 

- 

- 

NM 

DES J033520.97-540348.2 

53.83736 

-54.06338 

18.95 

18.46 

62.31 ± 1.06 

2.16± 0.26 

-2.50± 0.15 

1.15 ± 0.28 

-3.20 ± 0.30 

M 

DES J033524.00-540226.7 

53.85002 

-54.04075 

20.31 

19.89 

59.81 ± 1.83 

- 

- 

- 

- 

M 

DES J033531.14-540148.2 

53.87975 

-54.03007 

17.64 

17.07 

57.66 ± 0.96 

1.72 ± 0.20 

-3.02 ± 0.14 

1.88 ± 0.26 

-2.92 ± 0.16 

M 

DES J033533.71-535025.4 

53.89045 

-53.84038 

18.75 

18.25 

157.61 ± 2.22 

- 

- 

- 

- 

NM 

DES J033535.44-540254.9 

53.89766 

-54.04857 

20.68 

20.35 

61.74 ±4.79 

- 

- 

- 

- 

M 

DES J033536.94-535445.1 

53.90391 

-53.91253 

20.56 

20.16 

299.87 ±2.00 

- 

- 

- 

- 

NM 

DES J033537.06-540401.2 

53.90442 

-54.06701 

18.57 

18.03 

64.44 ± 1.14 

1.94 ± 0.21 

-2.70± 0.13 

1.67 ±0.48 

-2.86 ± 0.32 

M 

DES J033540.70-541005.1 

53.91957 

-54.16809 

18.23 

17.92 

45.83 ± 0.99 

- 

- 

- 

- 

NM 

DES J033544.18-540150.0 

53.93409 

-54.03056 

20.40 

19.91 

61.86 ± 1.99 

- 

- 

- 

- 

M? 

DES J033546.17-540733.9 

53.94235 

-54.12608 

19.23 

18.74 

216.54 ± 1.16 

3.51 ± 0.41 

- 

- 

- 

NM 

DES J033547.47-535926.5 

53.94781 

-53.99068 

21.25 

20.99 

-8.53 ± 4.78 

- 

- 

- 

- 

NM 

DES J033549.96-540321.5 

53.95816 

-54.05596 

19.66 

19.21 

289.73 ± 1.23 

4.93 ± 0.36 

- 

4.65 ± 0.59 

- 

NM 

DES J033550.10-540139.2 

53.95873 

-54.02756 

19.69 

19.26 

59.15 ± 8.25 

- 

- 

- 

- 

M? 

DES J033552.08-540733.9 

53.96700 

-54.12608 

20.56 

20.21 

152.93 ± 2.52 

- 

- 

- 

- 

NM 

DES J033556.28-540316.3 

53.98449 

-54.05452 

18.85 

18.37 

65.15 ± 1.18 

1.56 ± 0.24 

-2.88± 0.18 

1.24 ± 0.28 

-3.14 ± 0.27 

M 

DES J033558.15-540204.8 

53.99228 

-54.03466 

19.30 

18.82 

65.67± 1.12 

1.92 ± 0.38 

-2.57 ± 0.24 

1.33 ± 0.48 

-2.97 ±0.40 

M 

DES J033601.76-540405.5 

54.00734 

-54.06819 

19.56 

19.06 

63.25 ± 1.40 

1.08 ± 0.32 

-3.16± 0.37 

- 

- 

M 

DES J033606.25-540144.5 

54.02604 

-54.02903 

20.39 

19.93 

29.99 ±4.19 

- 

- 

- 

- 

NM 

DES J033607.75-540235.6 

54.03230 

-54.04321 

17.43 

16.82 

59.07 ± 0.92 

2.22 ± 0.21 

-2.78± 0.12 

- 

- 

M 

DES J033618.68-535745.1 

54.07784 

-53.96254 

18.04 

18.20 

55.71 ± 3.23 

- 

- 

- 

- 

M 

DES J033621.86-540040.7 

54.09109 

-54.01130 

20.30 

19.87 

64.66 ± 1.75 

- 

- 

- 

- 

M 

DES J033622.83-535955.5 

54.09513 

-53.99876 

18.58 

18.12 

126.61 ± 0.94 

5.25 ± 0.33 

- 

- 

- 

NM 

DES J033635.75-540156.9 

54.14895 

-54.03248 

18.05 

17.61 

131.78 ± 0.94 

- 

- 

- 

- 

NM 

DES J033635.78-540120.2 

54.14909 

-54.02228 

20.17 

19.73 

56.68 ± 1.91 

- 

- 

- 

- 

M? 






PTRAEEE _ 

















DES J033523.85-540407.5 

53.84938 

-54.06875 

16.45 

15.65 

65.69 

2.42 ± 0.19 

-2.90± 0.10 

- 

- 

M 

DES J033524.96-540230.7 

53.85400 

-54.04186 

19.69 

19.17 

67.99 

5.30 ± 0.57 

- 

- 

- 

NM 

DES J033539.04-535556.4 

53.91267 

-53.93234 

17.85 

17.16 

-20.65 

6.97 ± 0.26 

- 

- 

- 

NM 

DES J033548.04-540349.8 

53.95017 

-54.06384 

18.27 

17.69 

62.78 

3.45 ± 0.22 

-2.02 ± 0.11 

- 

- 

M? 

DES J033603.90-540026.5 

54.01625 

-54.00737 

17.79 

17.09 

-25.16 

5.87± 0.23 

- 

- 

- 

NM 

DES J033623.52-540518.9 

54.09798 

-54.08857 

19.75 

19.23 

213.51 

4.33 ± 0.55 

- 

- 

- 

NM 






_ PA/TOc; _ 











- Vj iVi w o - 






DES J033539.85-540458.1 

53.91605 

-54.08281 

18.33 

18.59 

69.30 

— 

— 

— 

— 

M 


^ Quoted magintudes represent the PSF magnitude derived from the DES coadded images by SExtractor. The one exception is DES J033523.85—540407.5 
which uses the median AUTO magnitude derived from the individual DES single epoch images. 
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Figure 6. (Left) Velocity as a function of radius for the M2FS spectroscopic sample. Ret II members are plotted as filled red symbols, 
and non-members are shown as small black dots. (Right) Metallicity as a function of radius for the GIRAFFE spectroscopic sample; 
symbols are the same as in the left panel. Note that CaT metallicity measurements depend on the absolute magnitude of the star, so for 
non-member stars the displayed points do not represent the true metallicity because those stars are not at the distance of Ret 11. 
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